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Abstract. We have used modern supercomputer facilities to carry out extensive 
Monte Carlo simulations of 2D hopping (at negligible Coulomb interaction) in 
conductors with the completely random distribution of localized sites in both space 
and energy, within a broad range of the applied electric field E and temperature T, 
both within and beyond the variable-range hopping region. The calculated properties 
include not only dc current and statistics of localized site occupation and hop lengths, 
but also the current fluctuation spectrum. Within the calculation accuracy, the model 
does not exhibit l/f noise, so that the low-frequency noise at low temperatures may 
be characterized by the Fano factor F. For sufficiently large samples, F scales with 
conductor length L as (L c /L) a , where a = 0.76 ± 0.08 < 1, and parameter L c is 
interpreted as the average percolation cluster length. At relatively low E, the electric 
field dependence of parameter L c is compatible with the law L c cx _E~ - 911 which 
follows from directed percolation theory arguments. 



PACS numbers: 72.20.Ee, 72.20.Ht, 73.50.Td 
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1. Introduction 

The theory of hopping transport in disordered conductors PJ El IS] at negligible Coulomb 
interaction is often perceived a well-established (if not completed) field, with recent 
research focused mostly on Coulomb effects. However, only relatively recently it was 
recognized that shot noise (see, e.g., Ref. [3]) is a very important characteristic of electron 
transport. In particular, the suppression of the current fluctuation density Si(f) at low 
frequencies, relative to its Schottky formula value 2e (I) (where (I) is the dc current), 
is a necessary condition for the so-called quasi-continuous ("sub-electron") charge 
transfer in such finite-current experiments as single-electron oscillations 

Earlier calculations of shot noise at hopping through very short samples (e.g., across 
thin films jZllHl) and simple lattice models of long conductors in ID [§| and 2D jTU] have 
shown that such suppression may, indeed, take place. However, calculations for the 
more realistic case of disordered 2D conductors have been limited to just one particular 
value of electric field ^Hj. It seemed important to examine whether the law governing 
this suppression is really as general as it seemed. Such examination, carried out in this 
work, has become practical only after the development of a new, advanced method of 
spectral density calculation ^T] in combination with the use of modern supercomputers. 
(The work reported below took close to a million processor- hours of CPU time.) 

As a useful by-product of this effort, we have obtained accurate quantitative 
characterization of not only the dependence of the average current on both temperature 
T and electric field E, but also the statistics of localized site occupation and hop lengths, 
which give a useful additional insight into the physics of hopping transport. 



We have studied hopping in 2D rectangular (L x W) samples with "open" boundary 
conditions on the interface with well-conducting electrodes ^Hj - see inset in Fig. [0 In 
the present study, we have concentrated on broad samples with width W 3> L c , where 
L c is the effective percolation cluster size (see below). The conductor is assumed to be 
"fully frustrated" : the localized sites are randomly distributed over the sample area, and 
the corresponding electron eigenenergies are randomly distributed over a sufficiently 
broad energy band, so that the 2D density of states uq is constant at all energies relevant 
for conduction. Electrons can hop from any site j to any other site k with the rate 



where r^ = y (xj — x k ) + (yj — y k ) is the site separation distance, and a is the 
localization radius ^2]- Such exponential dependence on the hop length has been 
assumed in virtually all theoretical studies of hopping. (The corrections to this law 
due to phase interference effects [T3 1 IT U ITK] are typically small.) However, in contrast 
to most other authors, we take Eq. literally even at small distances ~ a; this 
range is important only at very high fields and/or temperatures where the average value 



2. Model 




(1) 
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Figure 1. Linear conductivity a as a function of temperature T. Points show the 
results of averaging over 80 samples varying in size (L x W) from (20 x 10)a 2 to 
(160 x 100)a 2 , biased by a low electric field E <C Et- Points are Monte Carlo results. 
The thin dashed line is just a guide for the eye, while the thick solid line corresponds 
to the best fit of the data by Eq. ®. Here and below, the error bars are smaller than 
the point size, unless they are shown explicitly. (The error bars correspond to the 
uncertainty in averaging over an ensemble of independent samples that is larger than 
the calculational uncertainty for each of the samples.) The inset shows the system 
under analysis (schematically). 



of rjk becomes comparable to a. Of course our quantitative results for this particular 
region are only true for the localized states with exponential wavefunction decay. 

Another distinction from some other works in this field is that we assume that the 
hopping rate amplitude Tjk depends continuously on the localized site energy difference 
AU Jt = <f > - 4°> + e Er Jk : 

ff -W^ l- g pl-M iiW <2) 
This model coincides, for low phonon energies AUjk, with that described by Eqs. 
(4.2.17)-(4.2.19) of Ref. [2] for hopping in lightly-doped semiconductors, and of course 
satisfies the Gibbs detailed balance requirement Tjk = Tkj exp (AUjk/kBT). It is 
also close, but more physical than the "Metropolis" dependence which has a cusp at 
AU jk = 0. 

The interaction of hopping electrons is assumed negligible (with the exception 
of their implicit on-site interaction which forbids hopping into already occupied 
localized states). This assumption is well-justified for practically important materials, 
in particular very thin films of amorphous silicon, which is the major candidate 
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material for the implementation of sub-electron charge transfer components in single- 
electronic circuits Indeed, the relative strength of the Coulomb interaction may 
be characterized by a dimensionless parameter x = (e 2 //ta) x (z/o<2 2 ), where k is the 
relative dielectric constant J7j. For a film of thickness t ~ a, u may be estimated as 
tu, where v is the 3D density of states. For undoped amorphous Si, v is of the order of 
10 16 eV _1 cm~ 3 , and only special treatments may increase it to ~10 20 eV~ 1 cm -3 - see, 
e.g., Ref. JH]- As a result, for k ~ 10 and a ~ 3nm (both numbers are typical for the 
midgap states in Si), x is much less than unity for the entire range of v cited above, 
so that there is a broad range (A (lnT) , A (\nE) ~ 31n(x -1 )) of temperature T and 
electric field E where the Coulomb interaction is negligible [T7j . 

With this assumption, our model has only three energy scales: kgT, eEa, and 
(z/oa 2 ) . In other words, there are two characteristic values of electric field: 

E T = ^ and E = (3) 

ea ez/ cr 

We will be mostly interested in the case of low temperatures T < To, where 

1 

is the field-independent scale of temperature, so that the field scales are related as 
E T < Eq. (The only role of the dimensionless parameter g introduced by Eq. (j2J) is to 
give the scale of hopping conductivity a = g(e 2 /h). Coherent quantum effects leading 
to weak localization and metal-to-insulator transition are negligibly small, and hence 
the formulated hopping model is adequate, only if g -C 1.) 

The dynamic Monte Carlo calculations were carried out using the algorithm 
suggested by Bakhvalov et al. [T^j, which has become the de facto standard for 
the simulation of incoherent single-electron tunneling [20J. All calculated variables 
were averaged over the sample, and in most cases, over several (many) samples with 
independent random distributions of localized sites in space and energy, but with the 
same dimensionless parameters L/a, W/a, T/T , and E/E . We have used a new, 
advanced technique [TT] of the noise (current spectral density) calculation to save 
simulation time. The used supercomputer facilities are listed in the Acknowledgments 
section below. 



- t^- 2 (4) 



3. DC Current 



In order to understand the relation between our model and the prior results in this field, 
we have started from the calculation of dc current (/) as a function of T and E. If 
the electric field is sufficiently small (E <C Et), then the current is proportional to E, 
and the transport is completely characterized by the linear conductivity o = (I) /WE. 
Fig. ^ shows the calculated conductivity as a function of temperature T. In the region 
T ^ T this dependence follows the exponential T dependence of the 2D Mott law 
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\ 1/3' 



— « A (T, 0) exp 



(T, 0) ^ 



(5) 



where A (T, E) and -B (T, £7) are dimensionless, model-dependent slow functions of 
their arguments. We have found that our results may be well fit by Eq. (jSJ) with 
the following pre- exponential function: A (T, 0) = (23.4 ± 1.4) (T/T ) (0 - 68±0 - 04) , and 
constant B (T, 0) = 2.0 ±0.2 [21]. This latter result may be compared with the following 
values reported in the literature. In Ref. [22], B was analytically estimated to be close 
to 2.1. A different value, 3.45 ± 0.2 has been found by mapping a random 2D hopping 
problem to the problem of percolation in a system of linked spheres [2"3*1 12] . Finally, 
a close value 3.25 (with no uncertainty reported) has been obtained using numerical 
simulations of hopping on a periodic lattice, with a slightly different model for the 
function T(AU) [21] • The difference between our result and the two last values is 
probably due to the differences between details of the used models. 

At higher electric fields (E > Et), dc current starts to grow faster than the Ohm 
law, so that if we still keep the above definition of conductivity a, it starts to grow with 
E (Fig. |2I). At T -> 0, the results are well-described by the expression [231231231211123 



-»A(0, E) exp 
Co 



B(0,E)^ 



iV 1/3 ' 



(6) 



The data for not very high fields (Et «C E <C Eq) may be well fit by Eq. (JB*j) 
with constant B (0, E) = 0.65 ± 0.02 and pre-exponential function A (0, E) = 
(9.2±O.Q)(E/E ) {0 - 80±om) . (Note that the value B (0, E) = 1.27 given in Ref. [TO] 
corresponds to a different pre-exponential function used for fitting.) 

Finally, Fig.[2]shows that when the electric field becomes comparable with the value 
Eq defined by the second of Eqs. (JHJ), dc current and hence, conductivity, start to grow 
even faster than the exponential E dependence of Eq. (JBJ). 

4. Hopping Statistics 

In order to understand the physics of hopping in the three field regions better, it is 
useful to have a look at the statistics of localized site occupation and hopping length. 
We have found that for all studied values of E and T, the probability of site occupation 
closely follows the Fermi distribution with the local Fermi level 

/i(r) =hl- eEr (7) 
(where is the Fermi level of the source electrode) and some effective temperature 
T e g . Points in Fig. El show T cS as a function of electric field E for several values of 
physical temperature T. Dashed lines show the result of the best fitting of the naive 
single-particle master equation 



dt 



J d 2 r J de' exp (-1^ {-T(e-e'+ eEr) / (e) [1 - / (e')] 

+r(e'-e-eEr)/(e')[l-/(e)]} 
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Figure 2. Nonlinear conductivity a = (I) /WE as a function of electric field E for 
several values of temperature T. Each point represents data averaged over 80 samples 
of the same size, ranging from (20 x 14)a 2 to (1000 x 700)a 2 , depending on T and E. 
Points are Monte Carlo results. Thin dashed lines are only guides for the eye, while 
the thick solid line shows the best fit of the data by Eq. ©. 




Figure 3. The effective temperature T e g of current carriers as a function of electric 
field E for several values of temperature T. Closed points: Monte Carlo simulation 
results; dashed lines: master equation results. The solid curve marked T — is only a 
guide for the eye. 
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by a stationary Fermi distribution. Equation (jSJ) would follow from our model if electron 
correlation (in particular, percolation) effects were not substantial. In reality, we can 
expect the results following from Eq. (JHJ) to be valid only in certain limits [SO]- For 
example, in the low field limit with E — > 0, both methods give T e g = T. At higher fields 
the effective temperature grows with the applied field, which "overheats" the electrons. 
At very high fields (E/E > 0.3) both methods agree again and give 

k B T cS « CeEa, C = 0.71 ± 0.02. (9) 

(A similar result, but with C ~ 1.34, for our definition of a, was obtained by Marianer 
and Shklovskii .SH] for a rather different model with an exponential energy dependence of 
the density of states uq.) However, at intermediate fields typical of "high- field" variable- 
range hopping (Et -C E <C E ), the master equation still gives the same result (jOJ and 
hence fails to appreciate that in fact T e g is proportional to E 2 ^ 3 [32J. In order to explain 
this result, let us discuss the statistics of hop lengths (Figs. Eland El). 

The two- and one- dimensional histograms in Fig.0]show the probability density P of 
a hop between two sites separated by the vector Ar, and also the density Pa weighed by 
the factor — H^jl, where each H is the total number of hops (during a certain time 
interval) in the indicated direction, i.e jk = j —>■ k. The latter weighing emphasizes the 
site pairs (J, k) contributing substantially to the net hopping transport, in comparison 
with "blinking" pairs which exchange an electron many times before allowing it to 
advance along the field. It is clear that at relatively high temperatures or low fields 
(E <C Et) the non- weighed distribution should be symmetric (Fig. Hk). Figure EJl 
shows that in this case the one- dimensional probability density is well approximated 
by the Rayleigh distribution, P (r) oc r exp (— r/a c g), with a c g m a. However, the 
weighed hop distribution is strongly asymmetric even in the limit E — > (Fig. 0Jd). 
This asymmetry is even more evident at low temperatures or high fields (E 3> Et)', in 
this case the distribution has a sharp boundary (Fig. Hk). Figure EJl shows those cases 
where the ID histograms deviate substantially from the distribution predicted by the 
master equation - see the first of Eqs. (fTtjj). 

Figure El shows the r.m.s. non-weighed (r rms ) and direction-weighed (r rmds ) hop 
lengths, defined, respectively, as 

r ^ 2 = Y, j ,k r2 j k( H jk + H kj ) 

and 

2 _ Ej,fc r jfc \Hjk - H kj\ 

'raids — ^ I u it I \ / 

(that are of course just the averages of the histograms shown in Fig. HJ, as functions of 
applied electric field for several values of temperature. At T — > 0, hopping is strictly 
one-directional (i.e., if Hjk ^ 0, then H^j = 0), so that r rms and r rm d s are equal. In fact, 
simulation shows that in this limit both lengths coincide, at lower fields following the 
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iV 1/3 



r rms = r rmds = Da[ — ) , E T < E < £ , (12) 



with .D = 0.72 ± 0.01. (We are not aware of any prior results with which this value 
could be compared.) 

This scaling of r in the variable-range hopping region is essentially the reason for 
the scaling of T e fj mentioned above; in fact, the hopping electron gas "overheating" may 
be estimated by equating ksT c s to the energy gain eEr rm d s , possibly multiplied by a 
constant of the order of one. For the effective temperature, this estimate gives 



k B T cS = const x eEr rmds = GeaEl /3 E 2/3 = G ^^ ' , (13) 

u a 2 

in accordance with the result shown in Fig. El Our Monte Carlo simulations give 
G = 0.60 ± 0.02; we are not aware of any previous results with which this number 
could be compared. 

At higher fields (E > 0.1E ) the hop lengths start to decrease slower, approaching 
a few localization lengths a (Fig. Ek)- In this ( "ultra-high- field" ) region, the energy 
range ~ eEa for tunneling at distances of a few a is so high that there are always some 
accessible empty sites within this range, so that long hops, so dominant at variable-range 
hopping, do not contribute much into conduction. 

At finite temperatures, the most curious result is a non-monotonic dependence of 
the r.m.s. hopping length on the applied field - see Fig. At E — > 0, r rms has to be 
field- independent, and there is no scale for it besides a. (As evident as it may seem, this 
fact is sometimes missed in popular descriptions of hopping.) In order to make a crude 
estimate of r rms in this limit, one can use the master equation (jSJ). In this approach, 
at thermal equilibrium (i.e. at T independent of r), the hop length probability density 
P (r) can be found as 

P (r) = 2nr J de J de' exp f^j {T (e - e>) f (e) [1 - f (e')] 

+r(6'-6)/(e')[l-/(6)]} 

oc r exp ( — ) , (14) 
\ a J 

in a good agreement with the results shown in Fig. EJi for this case. From Eq. (fT4*|) . we 

get 

2X i/2 f / P^^dr 1 1 '~ 



,2/3 



V6a w 2.45a, (15) 



J P (r) dr 

in a good agreement with numerical data shown in Fig. 

A similar calculation for r rm( j s may be obtained by expanding the tunneling rate 
r(e — e' + eEr) in small electric field as T(e — e') + eErr'(e — e')- The result is 



P d (r) = eEr 2 J d(f>\cos4>\ J de J de' exp (—^j 
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x {V (e' - e) f (e') [1 - / (e)] + T' (e — e') / (e) [1 - / (e')]} , 



Pd (r) oc r exp I — 



— r 



^rmds 



J Pd (r) r 2 dr 



1/2 



12a « 3.46a. 



/ P d (r) dr 

The Monte Carlo data (Fig. |5J), however, differ from this result 
E -> 0, r rmds obeys the Mott law 



''rmds 



1/3 



+ /a, 



TcT ( 



o- 



(16) 

showing that at 
(17) 



with the best-fit values H = 0.52 ± 0.05 and / = 2.0 ± 0.1. In contrast, the function 
r rms (T) is rather far from Eq. (JUJ), because the Mott law refers to long hops responsible 
for transport (with the average approximately corresponding to r rm d s ), while r rms reflects 
the statistics of all hops. 



5. Shot Noise 



The current noise simulation has been limited to the case of zero temperature for two 
reasons. First, in the opposite limit (Et ^> E) current noise obeys the fluctuation- 
dissipation theorem; as a result, its low-frequency intensity can be found from cr(T), and 
hence does not provide any new information. Second, the calculation of the spectral 
density Sj (/) of current fluctuations with acceptable accuracy requires much larger 
statistical ensemble of random samples than that of the average current, which means it 
becomes increasingly difficult to extend it to finite temperatures even with the advanced 
averaging algorithm and substantial supercomputer resources used in this work. 

Figure El shows a typical dependence of the current noise spectral density Si, 
normalized to the Schottky value 2e (/), on the observation frequency uj = 2irf. One 
can see a crossover from a low-frequency plateau to another plateau at high frequencies. 
As the sample length grows, the crossover becomes extended, i.e. features a broad 
intermediate range oji < uj < u>h, just like in ID systems with next-site hopping [0J. The 
position of the high-frequency end u>h of this region can be estimated in the following 
way. 

In all single-electron tunneling systems, the high-frequency plateau is reached at 
frequency Uh close to the rate Y of the fastest electron hops affecting the total current 
j3 HOI EH ■ (For example, in systems described by the "orthodox" theory of single- 
electron tunneling, ujh ~ r max rs AU m3iX /e 2 R ~ 1/RC, where R and C are, respectively, 
resistance and capacitance of a single junction [331 ES]-) In our current case, this means 
that u h ~ (T jk ) max ~ ((} / fi) [At/jfc exp ( — fjfc/a)] max . For this estimate, AUji~ can be taken 
as fcsT e fj from Fig. El while according to the histograms shown in Fig. Efd), the length 
of shortest hops, still giving a noticeable contribution to the current, can be estimated 
as ~ r rm ds/2. For the case shown in Fig. |H1 (E/Eq = 8.75 x 10~ 3 ), these estimates yield 
AU jk /k B T k2x 10~ 2 , {r jk ) min /a w 2.5, giving finally u h /w ~ 1.5 x 10~ 3 , in a very 
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Figure 4. (a)-(c) Two-dimensional and (d) one-dimensional histograms of hop lengths 
for two typical cases: (a) and (b): T/T = 5 x 1CT 3 , E/E = 1.37 x 1CT 4 (E < E T ) 
and (c): T = 0, E/E = 1.09 x 1(T 3 (E T < E). The shade-coding in panels (a) and 
(c) corresponds to the probability P of hops with given Ar = (Ax, Ay), while that in 
panel (b), to the probability Pd weighed by factor \Hjk — Hkj\ - see the text. (Since at 
T = there are no backward hops, for the case shown in panel (c), P and Pd coincide.) 
Panel (d) shows P and Pd, averaged over all directions of vector Ar, for the low-field, 
intermediate, and high-field cases. Dashed lines show the distribution (| 1 fift given by 
master equation, for the best-fit values of parameter a e g. 

reasonable agreement with numerical results shown in Fig. |H1 (Note that this simple 
estimate, giving a length-independent value for Uh, is only valid for relatively long and 
broad samples.) 

At u — >• 0, a crossover to 1/f noise might be expected, because the discussion 
of this effect in some earlier publications EI] was apparently independent of the 
Coulomb interaction between hopping electrons. However, within the accuracy of our 
simulations, we could not find any trace of 1/ /-type noise for any parameters we have 
explored. This fact may not be very surprising, because all the discussions of the 1/f 
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Figure 5. R.M.S. hop length r rms (solid points) and the weighed average hop length 
fimds (open points) as functions: (a) of applied electric field E for several temperatures 
T, and (b) of temperature at E — > 0. Tilted straight lines in panels (a) and (b) show 
the best fits by Eqs. I|12|) and l|17|) . respectively, while the horizontal thin lines show 
the values following from the master equation. Curves are only guides for the eye. 



noise we are aware of require the presence of thermal fluctuations which are absent in 
our case (T = 0). 

Since the low-frequency spectral density is flat, at T = it may be considered as 
shot noise [SB] and characterized by the Fano factor |3j 

F== 2e(I) • (18) 
Similarly, in order to characterize the flat high-frequency spectral density, we may use 
the parameter 

= ft(/->oo) 

00 " 2e</> " 1 J 
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Figure 6. Spectral density Si of current fluctuations normalized to the Schottky noise 
2e (7) as a function of observation frequency u>, measured in units of u> = g/hiy a 2 , 
for several values of sample length L for T = and E/Eq = 8.75 x 1CP 3 . Small points 
show results for W/a = 60, while open squares are for W/a = 120 (at L/a = 60). 
Horizontal lines correspond to the Fano factor for hopping through short samples with 
one and two localized sites 0|H]. Curves are only guides for the eye. 

Figure El^a) shows the average Fano factors F and as a function of L for several 
values of the applied electric field, while Fig. |^b) shows that the same data can be 
collapsed on universal curves by the introduction of certain length scales: L c for F and 
L h for Foo. 

For the high frequency case, 

Foc= i-f) > L^L h , (20) 

where, within the accuracy of our calculations, (3 = 1. Such dependence could be 
expected, because the high-frequency noise at hopping can be interpreted as a result 
of the "capacitive division" of the discrete increments of externally-measured charge 
jumps resulting from single-electron hops through the system jH]. When applied to 
uniform (ordered) systems, these arguments always give the result F^ w l/iV^ with 
Nh = L/d being the number of electron hops (d the hopping length along the current 
flow) necessary to pass an electron through the system, regardless of hop rates jHElEl- 
For T — * in the case of disordered conductors, L h in Eq. (|2U|) may then be interpreted as 
the average hop length along current flow. This interpretation turns out to be correct. 
Indeed, Fig. |H] shows that the parameter Lh obtained from Eq. (|271jl scales with the 
electric field in a manner similar to x rm d s , especially at low fields, where it follows the 
variable-range hopping dependence of Eq. (fT2^1 . 
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Figure 7. Average Fano factor F and its high-frequency counterpart F^ as functions 
of sample length L normalized to: (a) the localization length a, and (b) the scaling 
lengths L c (for F) and Lt (for Foo) (see Fig. below), for several values of applied 
field at T = and W L c . Horizontal lines correspond to the average Fano factor 
for hopping via one and two localized sites p] |Hj • Straight lines are the best fits to the 
data, while dashed curves are only guides for the eye. 



The low-frequency value F, in the limit (L <C L c ), is weakly dependent on length 
and approaches F pa 0.7, which not surprisingly is consistent with the prior results for 
hopping via 1 intermediate site (F = 0.75) and 2 such sites [S] (F = 0.707). The 
results for long samples are much more interesting. We have found that they may be 
reasonably well fit with a universal dependence 

F = f^j , L » L c . (21) 

Here a is a numerical exponent; in the current study we could establish that 

a = 0.76 ±0.08. (22) 
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Equation ()21jl and the value for a are compatible with our previous results [TU] ct = 
0.85 ± 0.07 (for the same model, but just one particular value of E) and a = 0.85 ± 0.02 
(for nearest neighbor hopping on uniform slanted lattices). 

Figure |H1 shows the fitting parameter L c as a function of electric field E. In the 
variable-range hopping region, it may be fitted with the following law, 

( E \ M 

L c = Ja[-^\ , J = 0.04 ± 0.01, fi = 0.98 ± 0.08. (23) 



This law may be compared with the result of the following arguments. According to the 
arguments given in Ref. JUj, parameter L c may be interpreted as the average percolation 
cluster length (up to a constant of the order of 1). The theory of directed percolation 
jSll EH E] gives the following scaling: 

^{whsY- (24) 

Here (x) is the r.m.d.s. hop length along the field direction, while x c is its critical 
value. According to Ref. jH], the critical index 5\\ should be close to 1.73. Due to the 
exponential nature of the percolation, \(x) — x c \ ~ a, while (x) should follow a field 
scaling similar to that given by Eq. (|12|) . (Square points in Fig. |H] show that this is 
true for our simulation results as well.) Thus for sufficiently large (x) we arrive at 
Eq. (|2^j) with /i = | (l + 6\\) ~ 0.911. Equation (24) shows that this value is quite 
compatible with our numerical result, thus confirming the interpretation of L c as the 
average percolation cluster length. 

Note that in the variable-range hopping regime, L c has a different field dependence 
and is much larger than the average hop length. However, as the applied electric 
field approaches Eq, both lengths become comparable with each other and with the 
localization radius a. 



6. Discussion 



To summarize, our results for average conductivity and hop statistics are in agreement 
with the well-known semi-quantitative picture of hopping, including the usual variable- 
range hopping at low fields (E <C E?) and "high-field" variable-range hopping at 
E T <C E <C Eq. However, our supercomputer-based simulation has allowed, for the 
first time, a high-precision quantitative characterization of hopping, for a particular but 
very natural model. Moreover, our model also describes the "ultra-high field" region 
(E ~ Eq) where the variable-range hopping picture is no longer valid, since from most 
localized sites an electron can hop, with comparable probability, to several close sites. 
(In the last region, there are no clearly defined percolation clusters; rather, electrons 
follow a large number of interwoven trajectories.) 

Our simulations of shot noise at 2D hopping have confirmed our earlier hypothesis 
[TU] that in the absence of substantial Coulomb interaction, in sufficiently large samples 
(L, W ^> L c ) the Fano factor F scales approximately proportionally to 1/L - see 
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Figure 8. The values of Ly t and L c giving the best fitting of shot noise results by 
Eqs. (|20|l and (|21|l as functions of applied electric field E (open and solid circles, 
respectively). Open squares show the simple and direction- weighed average hop length 
along the applied field direction, defined similarly to Eqs. (|10|l and JUJ. Straight lines 
are the best fits to the data. 

Eq. (j2"Tj) . Other confirmations of this hypothesis have come from recent experiments 
with lateral transport in SiGe quantum wells j32] and GaAs MESFET channels (321 El- 
Unfortunately, these experiments are not precise enough to distinguish the small 
difference between the exponent a in Eq. (j21j) and unity. 

The hypothesis that a is in fact equal to 1 for sufficiently long samples seems 
appealing, because it would mean the simple addition of mutually-independent noise 
voltages generated by sample sections connected in series. On the contrary, a deviation 
of a from unity would mean that some dynamic correlations of electron motion persist 
even at L ^> L c . For ID hopping this fact is well established: in at least one exactly 
solvable model (the "asymmetric single exclusion process", or ASEP j3S|) the dynamic 
correlations may change a from 1 to 0.5. However, for 2D conductors the fact that 
the correlation length (if any) may be substantially larger than the percolation cluster 
length comes as a surprise [46J. 

From the point of view of possible applications in single-electron devices ^Hj, the 
fact that F may be suppressed to values much less than unity is generally encouraging, 
since it enables the use of circuit components with quasi-continuous charge transport. 
However, in order to achieve the high-quality quasi- continuous transfer (say, F < 0.1), 
the sample length L has to be at least an order of magnitude longer than the percolation 
cluster length scale L c . On the other hand, L c itself, especially in the most interesting 
case of low applied fields, is substantially longer than the localization length a (Fig. (SJ), 
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which is of the order of a few nanometers in most prospective materials, e.g., amorphous 
silicon. Hence, it may be hard to implement sub-electron transport in conductors 
substantially shorter than ~100 nm. This size is quite acceptable for experiments at 
low temperatures (say, below 1 K), but is too large for the most important case of 
room-temperature single-electron circuits ^E] , because of large stray capacitance of the 
hopping conductor, which effectively sums up with the capacitance of the island to be 
serviced. 
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